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Abstract 

X-ray absorption and its dependence on the polarization of light is a powerful tool 
to investigate the orbital and spin moments of magnetic materials and their orien- 
tation relative to crystalline axes. Here, we present a program for the calculation 
of dichroic spectra from first principles. We have implemented the calculation of 
x-ray absorption spectra for left and right circularly polarized light into the Wien2k 
code. In this package, spin-density functional theory is applied in an all-electron 
scheme that allows to describe both core and valence electrons on the same footing. 
The matrix elements, which define the dependence of the photo absorption cross 
section on the polarization of light and on the sample magnetization, are computed 
within the dipole approximation. Results are presented for the £2,3 and M4 5 egdes 
of CeFe2 and compared to experiments. 
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1 Introduction 



Dichroism is the property of a material to absorb photon beams of different 
polarization with different cross-sections. In particular, x-ray magnetic circular 
dichroism (XMCD) originates from excitations of core electrons to unoccupied 
spin-split conduction states, giving rise to different absorption behavior of left 
and right polarized light. XMCD, together with x-ray absorption spectroscopy 
(XAS), is one of the most effective tools for obtaining information about mag- 
netic systems. It exhibts different capabilities of characterizing these systems 
which cannot be afforded by traditional magnetic techniques. Among the ad- 
vantages of XMCD, it is worth mentioning its high sensitivity and chemical 
selectivity, the latter being an essential property to study magnetism in alloys, 
oxides, impurities, surfaces, and interfaces. Moreover, it represents one of the 
few methods which can discriminate between orbital and spin contribution 
to the total magnetic moment. Information about spin and orbital magnetic 
moments, associated to the photo-absorbing atom are extracted from XMCD 
spectra, exploiting two specific sum rules [Tf2] which allow to extract the ex- 
pectation value of the orbital and spin angular momentum operators, L z and 



We have developed a package for the calculation of x-ray absorption spectra 
at the K-, L-, and M-edges for left and right circularly polarized light in the 
Wien2k code [3] More specific, we have extended the Wien2k OPTIC package 
[3], which allows to calculate optical properties. The matrix elements, which 
define the dependence of the photo-absorption cross-section on the polariza- 
tion of light and on the sample magnetization, are calculated within the dipole 
approximation. 

In the following, we will give an extensive description of the formulae be- 
hind the method. As prototypical examples, we have investigated dichroic 
and absorption spectra at the M 4 5 and L 2) 3 edges of cerium in CeFe^ which 
demonstrate the applicability of our code to such core excitations. 



2 XAS calculations can already be performed in the Wien2k distribution via the 
XSPEC package, but in that implementation the polarization of light is not con- 
sidered, thus dichroic spectra can not be trivially derived. 
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2 Theoretical background 

XMCD is determined by the difference in absorption between left and right 
circularly polarized light, where left and right are referred to the propagation 
direction of incoming radiation with respect to the magnetization axis of the 
system. In particular, photons are left (right) polarized or, equivalently, their 
helicity is — H when the direction of propagation is anti-parallel (parallel) 
with respect to the magnetization. 

The absorption of a photon with polarization vector e results in the excitation 
of a selected atom from an initial state into a final state |^/). In the 
framework of the electric-dipole approximation, the cross section fi(u) for x- 
ray absorption is given by Fermi's golden rule: 

f i(u)<x y £ l \e--D fi \ 2 6{E f -E i -hu)), (1) 
/ 

where 

D/i =< tt/lpltti > (2) 

is the matrix element of the momentum operator p = — ihV r between the 
initial and final states of energy Ei and Ej. It is convenient to express the 
operator V in terms of spherical components 

Similarly, the polarization vector for right (+1) and left (—1) polarized light 
is 

e±i = e x ± ie y (4) 

and the scalar product appearing in Eq. [IJwill select the ±1 components of 
the vector D/j. The absorption cross section for the two polarizations then 
becomes 

^(u) oc J2 \Dfl\ 2 S(E f -Ei- hw), (5) 
/ 

where 

^ 1 = 6 ± -(* i |p|*/) = (* i |V±|* / >. (6) 

2.1 Sum rules 

A connection between the integrated XMCD signal and the ground state ex- 
pectation value of the projection of orbital angular momentum on the magne- 
tization axis was suggested by Thole, Carra, and van der Laan pp. Later, the 
same authors derived a second sum rule for the spin moment using graphical 
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angular momentum techniques [2]. They developed these sum rules by ana- 
lyzing near-edge x-ray circular dichroism and using a single-ion model with 
a partially filled valence shell. Starting from the absorption cross section for 
pure dipole transitions, Eq. (TJ integrating over the photon energy and applying 
the Wigner-Racah coupling techniques [5f6]T7f8] . they obtained the following 
relations: 

I 3± du(fi + -^) i c ( c+ i)-l(l+i)_2 
J j± dw{n+ + /I- + ijP) 2 £{£ + l){A£ + 2-n) [ z) U 



Sj+ duj (^ + ~ V ) - [(c+ l)c] - /i ) 



-2-c(c+l) 



+ 



3c(4£ + 2 - n) 
£(£ + l)[£(£ + l) + 2c(c 



(&>+ 
l)+4] 



3(c-l) 2 (c + 2)' 



6£c(£ + l)(4£ + 2-n) 



where c represents the core-hole orbital quantum number, n the number of 
electrons in the valence shell, I the valence orbital quantum number, j± = 
c ±1/2 the quantum number of the two partners of the spin-orbit split inner 
shell, fi = (/i + +/i_)/2 the absorption cross section for incident light polarized 
along the direction of magnetization, and T z the magnetic dipole operator 
defined as: 



T 7 , 



3r, : 



(9) 



3 X-ray circular dichroism within the LAPW basis set 



3.1 Density functional theory and the LAPW basis set 



In DFT-based band structure calculations [9], one-electron wavef unctions, 
\l/ n (r), and eigenvalues, E n , are derived solving the Kohn-Sham(KS) equa- 
tions fTOl: 



-V 2 + V H {r) + V nucl {*) 



+ V xc (r) 



tf n (r) = EA(r) 



(10) 



where V#(r) denotes the Hartree potential, V nuc i{r) is the bare Coulomb po- 
tential of the atomic nuclei, and V xc (r) is the exchange-correlation potential, 
which is defined as the functional derivative of the energy E xc [n], n(r) being 
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the ground-state density of the system. In order to solve these equations nu- 
merically, the KS orbitals are expanded in terms of an appropriate finite set 
of basis functions {4> u }, 

*n.(r) = £C<>,(r), (11) 

v 

and the coefficients C™ are then obtained by diagonalizing the Hamiltonian 
matrix. 

If one is interested in exploring the physics of the core region, as in the case 
of the XMCD, all-electron schemes are needed. Among them, the ones de- 
scending from the augmented planewave (APW) method are the most precise 
ones for a proper description of magnetic properties. These methods employ 
a hybrid set of basis functions, i.e., atomic-like basis functions in the unit-cell 
region close to nuclei and planewaves elsewhere. Our method of choice is the 
linearized augmented planewave (LAPW) method [TTlll2lll3|rH] . where the ba- 
sis functions inside the muffin-tin spheres are linear combinations of the radial 
functions Ug(r,Eg) and their energy derivatives ill (r, Eg) at the trial energy 
Ei times spherical harmonics Yg m (r). 

^ G (S a |r)^[^(k+GKM f )+B;(k+G)^M# k (f) (12) 

£m 

with G denoting a reciprocal space vector. The radial functions are obtained 
by solving the radial Schrodinger equation in the spherical potential of the 
respective atomic sphere. The coefficients A" m (k + G) and Bg m (k + G) are 
determined for each atom by matching the two types of basis functions at the 
atomic sphere boundary, both in value and slope. Fixing the energies Eg in 
Eq. [12] makes the basis set energy independent. As a consequence, the secular 
equation becomes linear in energy, leading to a generalized eigenvalue problem. 
But, at the same time, one faces the restriction that for a given I value only 
the sates of one principal quantum number can be described. A solution of 
this issue is the introduction of additional basis functions, called Local Orbitals 
PS] (LO) of the form 

0lo(S q + r) = [A^uftr, Eg) + B&,u?(r, Eg) + Q>°(r, E lo )] Yg m {v) (13) 

where Eg is the same as in the LAPW basis (Eq. [T2"|) . and E\ Q represents the 
trial energy of the semicore state. Here A^ m and Bf m are determined such 
that the LO and its slope go to zero at the sphere boundary, i.e., these basis 
functions are completely confined within the atomic spheres, whereas Cf m are 
chosen in order to normalize the basis function. 

An alternative to the linearization described above is provided by the APW+lo 
method [IB]. Here, the corresponding basis consists of APW functions taken 
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at a fixed energy Ei 

k+G (S a + r) =J^(k + GK(r,£i)y tol (f) (14) 

supplemented by local orbitals of the form 

MS* + r) = [i^°<(r, JSSe) + ^°<(r, Y im (v). (15) 

The two coefficients and are determined by normalization and by 
requiring local orbitals to have zero value at the muffin-tin boundary, but not 
zero slope. 

The Wien2k code can supply both types of basis functions, i.e., APW+lo and 
LAPW+LO, and it is possible to choose one or the other independently for 
different atoms and angular momentum numbers. 

3.2 The XMCD formalism in the LAPW basis 



As we have seen in Sec. [2J the main ingredient for the calculation of a dichroic 
signal is the expression of the x-ray absorption yield. Within the one-particle 
framework and in the dipole approximation, the absorption cross-section /i for 
incident x-rays is determined by the probability of an electron to be excited 
from a core state ^^ mj {y) with energy Ej, to a final valence state ^(r) with 
energy E^, according to Fermi's golden rule (Eq. Note that core states are 
characterized by atomic quantum numbers j, rrij and valence states by band 
index / and momentum k. From this formula, dichroic and total absorption 
signals can be calculated through the linear combinations fi + —fi~ and fi + +fi~ , 
respectively. Therefore, the main task is to compute the momentum matrix 
elements. 

In Wien2k, core states are calculated by the routine LCORE, which represents 
a modified version of the relativistic atomic LSDA code by Desclaux [TTffTS] 
that solves the fully relativistic Dirac equation. Core states wave function can 
thus be written as: 

Wr)= E E ^£ tSm y ML (v), (16) 

lf=-I» J= -l/2 

where 

= <(r)Y LM (r) (17) 

and u c L (r) are the solutions of the radial part of the Dirac equation. Valence 
states are expanded in terms of the LAPW basis set according to Eq. [UJ here 
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rewritten as: 

v£/ ?W = EC , / k (G)0 k+G (r) 

G 

Inserting Eqs. [16] and [18] in Eq. [6] we obtain 



EEC , / k (G)(E^ ms ^MW|V ± |0 k+G (r)}. (19) 



G M 



These matrix elements must be inserted into Eq. [5] which, after summing over 
rrij, becomes: 

E \V%, mj (V\ 2 5(E k f -E 3 -hto) (20) 
The V± operators can be expressed in spherical coordinates as 

V ± = sin9e^ + -e^ (cos9^- ± -4^) ■ (21) 
or r \ m svnU o<p J 

In the case of positive helicity of the photon, the V+ operator applied to the 
valence state gives: 

(V + )0 k+G (r) = 

= E ^ Kn(k + GK(r, E t ) + Bf m (k + G)uf(r, E t )] sin6e^Y im (r) + 

tva 

+ \ E + GK(r, E t ) + B? m (k + G)«?(r, E e )] ■ 



(22) 



Exploiting the following relations between spherical harmonics: 

e^sinBY lm = F^Y £+1>m+1 + F®Y £ _ 1>m+1 (23) 

e- i(t >sineY lm = F^Y t+ltm . 1 + F^Y e ^, m -i (24) 

^ ( cos6 §q +^0^l) Y ^ = -^iS*Ww + (' + 1)^2^-1^+1 (25) 



with 



r tm 



(l + m + l)(l + m + 2) 
\ (2£ + l)(2£ + 3) 1 j 



(e-m)(e- m -i) 

\ (2£ - 1)(2£+ 1) 1 ' 
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(£-m + l){£-m + 2) 
\ (2£+ l)(2£ + 3) 



(2g 



P (4) 



(£ + m)(£ + m - 1) 



(29) 



^ (2£-l)(2l+l) 
the momentum matrix element relative to right polarized radiation becomes: 

D^, mj (k) = (^, mj |V + |^) = 

= EEE C /( G ) / r 2 dr [ dn^ClZs m A(r)Y£ M (ry 



G M Im 



^{[AiimU'z + B^ m M^ -- (A im Ui + Si m ^)]F^l£ + i jm+1 + 



m« ' rs " 



+ [A £m u' £ + i?£ m M£ ^ (A £m w^ + Bi m ut)]F^Ye-i, m +i} 



r v- 



EEEE^,«n.Cf(G){ / ^HIC - ^]F« / dUY* LM Y t+1 



m-t-1 



+ 



G M £m rn s 

+ J r 2 dru c L [WL + Ltlwt^F® J dnY£ M Y £ _ lj7n+1 } = 

v * ' 

L 1/2 

= E^/ k (G) E E c£2U 

G M=-Lm s =-l/2 

{A i _ 1)M _ 1 [ y uK-ir 2 * - (L - 1) y wiw i _irdr]F^ liM _ 1 + 
+ B L _ 1M -i[J u c L u' L ^r 2 dr — (L — 1) y u£u i _irdr]F| 1 2 1;M _ 1 + 
+ y nK +1 r 2 cir + {L + 2) J u c L u L+l rdr}F[ 2 l hM _ 1 + 

+ B L+1>M ^[J u c L u' L+1 r 2 dr + (L + 2) J ulu^rd^F^^} 



SM,m+l 



(30) 



The corresponding momentum matrix element related to the V_ operator is 
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derived in the same way, and the following formula is obtained: 



L 1/2 



% m ,(k) = E*) E E c i 



G M=-L m s =-l/2 



{A L _ M/ +i[ | ^<_ x r 2 dr ~{L-l)J u c L u L ^rdr]Fi 3 \ M+1 + 
+ B L _ ltM+1 [j u c L v! L _ x r 2 dr - (L - 1) | wi«L-irdr]Ffi ljM+1 + 
+ A L+lM+1 [J u c L u' L+1 r 2 dr + (L + 2) J u c L u L+1 rdr]F^l lM+l + 
+ #z+i,m+i[ / <M' L+1 r 2 rfr + (L + 2) | w^irdr]^^} 



(31) 



Finally, the momentum matrix elements, Eqs. [3U] and EH have to be inserted 
in Eq. [201 an d the absorption spectra are then calculated for the two polariza- 
tions by summing over band indices and momenta. The total absorption and 
dichroic signals are derived as sum and difference of fi + and respectively. 

A few more words should be spent to see what happens when Local Orbitals 
are taken into account. It is worth noticing that, in this case, the basis set is 
augmented by a term given by EqJTBl The Kohn-Sham wavefunction (Eq. [TBI) 
thus becomes 



This means that a further term must be added to the expression for the 
matrix elements (Eqs. 13111 and I3"T|) . As, however, the derivation is along the 
same lines as above, it will not be explicitely described. The same reasoning 
can be applied, mutatis mutandis, if the APW+lo basis set (Eq. [15]) is used. 

If we assume the core-hole and the excited photoelectron to have both an infi- 
nite lifetime, the initial and final energies of the electronic transition are sharp, 
resulting in rapidly varying spectra. In reality, however, initial and final states 
have finite lifetimes: The core-hole decays by radiative or Auger electronic 
transitions, whereas the excited electron can lose energy by emitting plas- 
mons or creating electron-hole pairs. These finite lifetimes can be accounted 
for by a smearing of the spectra, i.e., adopting a Lorentzian broadening for the 
initial and final states. Usually, the (smaller) energy-dependent photoexcited 
electron broadening can be safely neglected, and only the more effective core- 
hole broadening (r c ) is considered. The broadened spectra, F(E), is obtained 
by F(E), the unbroadened one, by: 



Particular care must be taken in choosing the broadening parameters. In the 
case of the L 2 ,3 edge, for instance, different values must be used for the two 




(32) 




(33) 
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edges because, in the experiments, the L2 line is wider than L3. This different 
broadening is determined by the fact that the decay of a px/ 2 core- hole has a 
channel (the super- Coster-Kronig process pi/2 — > P3/2) which is not available 
to a P3/2 core-hole. Finally, the spectra can be convoluted with a Gaussian 
function in order to take into account the finite spectrometer resolution. 



4 Implementation into the Wien2k code 

The application of our package is based on a self-consistent calculation for a 
spin-polarized system, including spin-orbit coupling. It concerns two programs, 
which are OPTIC and JOINT, that must be run consequently. The former per- 
forms the calculation of the momentum matrix elements m .(k) (Eq. [T9|) . 
whereas JOINT carries out the integration over the Brillouin zone (BZ) by 
means of Bloechl's tetrahedron method |19j . 

OPTIC reads the core and valence wavefunction, which are calculated by the 
routines LCORE and LAPW1, from the files case . vectorsoup/dn and case . corewf up/dn, 
respectively. In order to perform the XMCD calculation, the user has to de- 
fine the atom for which the spectra should be calculated as well as the kind 
of transition. These parameters must be provided in the input file case . inop, 
which is described in Table [TJ 

Table 1 

A typical input file case . inop for computing the matrix elements which determine 
the dichroic spectra. The values correspond to the calculation of the Ce 1^2,3 edge, 
shown in Sec. 15.11 

99999 1 KUPLIMIT, KFIRST 

0.0 3.2 EMIN, EMAX 

1 1 L23 XMCD, ATOMJJUM, EDGE 

The meaning of the input parameters is the following: 

• KUPLIMIT is the maximum number of k-points to be taken into account; 
KFIRST is the first k-point. (This value should be only different from 1 for 
special analysis purposes.) 

• EMIN and EMAX define the absolute energy range (in Ry) in which the matrix 
elements are calculated. 

• XMCD is a parameter to decide whether to perform an XMCD calculation 
(XMCD = 1) or optical properties with the original version of the OPTIC 
program (XMCD = 0). 

• AT0M_NUM specifies the atom (according to case . struct file) for which an 
XMCD calculation is performed. 
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• EDGE indicates the particular edge; possible values are (case sensitive): Kl, 
LI, L23, Ml, M23, and M45. 

The most important routines of the OPTIC program, which are relavant for 
the XMCD calculation are listed below together with a short description: 

• opmain.f is the main routine, reads the input parameters from the file 
case . inop and the structural data from case . struct and calls sph-UPcor . f rc. 

• sph-UPcor . frc is the XMCD master routine. It is responsible for the com- 
putation of the momentum matrix elements according to Eq. [191 

• momradintc . f calculates the radial integrals that appear in Eqs. 1301 and I3T1 

• mmatcv.f performs the summation of the radial integrals over M (Eqs. [301 
and EH). 

• outmat ABZ . f calculates the squared matrix elements and sums over the the 
quantum number hl, (Eq. [201) . 

The momentum matrix elements between the selected core state and the con- 
duction states are stored in case . symmatlup for the upper spin-orbit split 
core state (e.g. L 3 ) and in case . symmat2up for the lower spinorbit-split core 
state (e.g. L2) for each k-point (k) and every band (f). In the case of K, Li, 
and Mi edges, only case . symmatlup is writtenjjO 

As mentioned above, the integration over the BZ is carried out by JOINT with 
the input file (case . injoint) described in Tabled 

Table 2 

A typical input file case . injoint for computing the dichroic spectra. The values 
are the ones used for the calculation of the Ce L2 3 edge, shown in Sec. 15.11 



1 9999 9999 NYMIN, NYMAX, NYOCC 

-1.0 0.001 2.0 EMIN, DE, EMAX 

ryd ECV 

1 XMCD 

-412.2 -445.0 C0RE_ENE1, C0RE_ENE2 

3.4 3.2 C0RE1BRD, C0RE2BRD 

0.1 SPECTR 



The meaning of the input parameters is the following: 

3 Core states are treated fully relativistically. As a consequence, those states with 
principal quantum number £ different from zero are spin-orbit split, resulting in a 
state whose total angular momentum J is given by the sum of orbital and spin 
momentum (J = I + s), and in another state in which I and s couple antiparallel 
(J = I — s). As an example, a 3d core state (£ = 2) will split in 3d (J = 5/2) and 3d* 
(J = 3/2) states, with the energy of the former higher than the one of the latter. 
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• NYMIN, NYMAX are the lower and upper band index for the summation (Eq. 
[20]) . Optionally, one can provide the index of the highest occupied band 
NYOCC (known from the ground-state calculation) to avoid summation over 
terms (of zero value) between occupied states. 

• EMIN and EMAX define the energy window (in Ry) with respect to the Fermi 
level, for which the spectra should be printed; DE is the corresponding in- 
crement for this energy mesh. 

• ECV defines the units for the output (case sensitive): 'eV, 'ryd', or 'cm-l' 
for eV, Ry, and wavenumbers (cm -1 ), respectively. 

• XMCD is the parameter to distinguish between an XMCD calculation (XMCD 
= 1) or a calculation of optical properties (XMCD = 0). 

• C0RE_ENE1 and C0RE_ENE2 are the energies of the higher (Is, 2s, 2p, 3s, 3p, 
or 3d) and lower (2p*, 3p*, or 3d*) core states, respectively. In the case of 
Ki, Li, and M 1 , C0RE_ENE2 is not read. 

• C0RE1BRD and C0RE2BRD are the core-hole lifetime broadenings in eV (i.e., 
T c in Eq. [33]) of C0RE_ENE1 and C0RE_ENE2, respectively. Again, in case of 
K\, Li, or Mi core2brd is not read. 

• SPECTR is the spectrometer broadening. 

The momentum matrix elements are read by the JOINT routine readopx.f 
from the files case . symmatlup and/or case . symmat2up (depending on the 
edge), and then the BZ integration is performed by the main routine joint . f 
and the subroutine arbdosx.f. Broadened and unbroadened dichroic and ab- 
sorption spectra are then written into case.xmcd and case .rawxmcd, respec- 
tively. To give the user the choice to apply sum rules (see Sec. 12. ip . broad- 
ened and unbroadened (raw) spectra for the single edges are written to the 
files case .broad2p, case .broad2ps, case.raw2p, and case.raw2ps, where 
2p refers to the upper edge and 2ps to the lower one. 



5 Results 



XMCD spectra can be calculated for any kind of system for which a Wien2k 
ground-state calculation has been performed, such as bulk materials or su- 
percell calculations. In Refs. [2"Uf2"Tj results for bulk iron, cobalt, and nickel 
are reported and compared to experimental as well as theoretical results from 
the literature. Overall good agreement was found in all cases. Here, we show 
two case studies: an example of a transition from c?-like core states to /-like 
valence bands, in particular the edge of cerium in CeFe2, and a dichroic 
excitation at the L 2) 3 cerium edge of the same system. 
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5.1 CeFe 2 



This system has been extensively investigated both theoretically [23|24|25] 
and experimentally [26|27|28||29] . where the main interest was focused on its 
anomalous physical properties in comparison with other rare-earth (R) RFe2 
compounds. In particular, it exhibits very low magnetic moments at iron and 
cerium sites (with experimental values of 1.40 \ib and 0.70 respectively) 
which are antiferromagnetically coupled to each other. 

CeFe2 crystallizes in the cubic Fd3m structure with a lattice constant a=7.3lA. 
The band structure calculation has been performed by sampling the whole 
Brillouin zone with 2000 k-points (163 in the irreducible part). Exchange and 
correlation effects are accounted for by the generalized gradient local approx- 
imation (GGA) in the version proposed by Perdew et. al. |22j . Spin-orbit 
coupling in the valence shell is included via a second variational scheme, as 
implemented in Wien2k. 

In Fig. [1], absorption and dichroic spectra are plotted for the edge of 
cerium. The theoretical absorption spectra (full line) present a very simple line 
shape at both the M 5 (826 eV) and M 4 (845 eV) edges, without any feature, 
whereas the experimental ones (dotted line) show a characteristic second peak 
at about 5 eV above the principal ones. This satellite has been discussed in 
the literature [28] and can be reproduced when including multiplet effects 
in the calculation. We will not focus on this point here, since it is outside 
the scope of this work. Lorentzian broadenings of 1.4 eV and 0.6 eV have 
been chosen for the M 5 and M 4 edges, respectively, in order to reproduce 
the experimental branching ratio at the two absorption edges. This choice 
leads to a good agreement also for the XMCD spectrum (right panel of the 
figure), apart from a small negative peak in the experimental M5 edge at about 
830 eV, related to the satellite, which is not reproduced in this calculation. 
Applying the sum rules, we have calculated the contribution of the / states 




820 830 840 850 860 820 830 840 850 860 
Energy [eV] Energy [eV] 

Fig. 1. Absorption (left) and dichroic (right) spectra for the M4 5 edge of cerium. 
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to the total spin and orbital moment to be [i SV in= -0.46 [Ib and /i rfc= -0.03 
Hb, respectively. 

As an example of a transition from p core states to d valence bands, we report 
in Fig. [2]Ce L2 and L3 XMCD spectra. Again, they are in very good agreement 
with the experimental ones. In this case, the two edges have to be depicted 
in two separate windows, because the spin-orbit separation between the two 
core states is very strong, i.e., about 32.8 Ry. 




Fig. 2. XMCD spectra of the Li (left) and the L3 (right) edge of chromium. The 
L>2 (L3) edge has been convoluted with a 3.2 eV (3.4 eV) full width half maximum 
Lorentzian. 



Applying sum rules to the L 2 3 spectra we found a contribution from the cerium 
d states to the total spin and orbital moment of cerium of -0.03 \ib and 0.00 
fis, respectively. We obtain a total magnetic moment of -0.52 \xb at the cerium 
site. This value is slightly smaller than the one (-0.69 \ib) found by Antonov et 
al. [30] with an atomic sphere approximation fully relativistic linear muffin-tin 
orbital (ASA-LMTO) code. In particular, they reported almost equal contri- 
butions from Ce / and Ce d to the total magnetic moment of Ce, whereas 
in the present work the Ce-/ contribution dominates the total moment. This 
small discrepancy can be explained by the different band-structure method 
and exchange-correlation potential adopted. As the expectation values of the 
moments depend in both methods on the atomic-sphere radii used, a slight 
difference may arise from this fact. 

In order to compute the total magnetic moment of the system, we have applied 
sum rules to the L 2i 3 spectra of iron (not shown here), finding 0.04 \xb and 1.25 
Hb for orbital and spin moments, respectively. By summing up all moments at 
the iron and chromium sites (fi tot =fi orb (Ce)+^ spin (Ce)+2[fi orb (Fe)+fi spin (Fe)}), 
we obtain a value of 2.08 /xg, which is about 9% smaller than the experimental 
value of 2.30 hb- 
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Another point worth to mention is that from the inferred moments of Fe 
(1.29 hb) and Ce (-0.52 fi B ), this calculation correctly predicts an antiparallel 
alignment between Fe and Ce, the same as found in experiments |28] as well 
as in earlier calculations |23]. 



6 Conclusions 

We have introduced a tool for the calculation of x-ray absorption spectra at 
the K-, L-, and M-edges for left and right circularly polarized light in Wien2k 
code, by generalizing the already implemented package OPTIC. As an example 
we have shown XMCD and XAS spectra of CeFe2 at the cerium M 4 5 and L 2; 3 
absorption edges. The results exhibit good agreement between theoretical and 
experimental spectra and demonstrate the capability of the code to reliably 
treat such core-level excitations. 
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